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ABSTRACT 






Assume a technology that permits undistorted laser beam 
propagation from the aft section of a streamlined turret. 

A comparison of power on a distant airborne target is made 
between a single aperture in a large scale streamlined turret 
with a turbulent boundary layer and various arrays of aper- 
tures in small scale streamlined turrets with laminar flow. 
The array performance is mainly limited by the size of each 
aperture. From an array one might expect, at best, about 
40 percent as much power on the target as from a single 
aperture with equal area. Since the turbulent boundary layer 
on the large single- turret has negligible effect on beam 
quality, the array would be preferred (if all development 
efforts were esentially equal) only if a laminar wake is an 
operational requirement. 



4 



TABLE OF CONTENTS 



Page 

I. INTRODUCTION 12 

II. METHOD-- - 15 

III. DIFFRACTION OF LASER BEAMS 17 

A. SINGLE APERTURE 17 

B. ARRAY OF APERTURES - 17 

C. IRRADIANCE ON A TARGET FROM AN ARRAY OF 

APERTURES 19 

IV. FLUID DYNAMICS RELATED TO TURRETS 30 

A. SEPARATED FLOW REGIONS 50 

B. WAKES FROM TURBULENT BOUNDARY LAYERS 3 0 

C. TRANSITION TO TURBULENCE IN THE BOUNDARY 

LAYER 3 0 

V. INFLUENCE OF TURBULENT BOUNDARY LAYER ON POWER IN 

THE BUCKET 54 

A. SELECTION OF MODEL 34 

B. ESTIMATE OF PIB FOR SAT 34 

VI. CONCLUSIONS - 39 

APPENDIX A PROGRAM LISTING GRID SUMMATION 41 

APPENDIX B PROGRAM LISTING TRAPEZOIDAL RULE 45 

LIST OF REFERENCES 45 



DISTRIBUTION LIST 



47 



LIST OF TABLES 



Table Page 

1. KINEMATIC VISCOSITY VS. ALTITUDE IN THE 

ATMOSPHERE 31 

2. FPIB FOR SELECTED VALUES OF D s /a AND aL 35 



6 



LIST OF FIGURES 

Figure Page 

1. PRESENT AND FUTURE TECHNOLOGY TURRET DESIGNS 14 

2. COORDINATE SYSTEM FOR AN ARRAY DIFFRACTION PROBLEM--- 18 

5. FPIB FOR VARIOUS ARRAY CONFIGURATIONS 23 

4. IRRADIANCE ON THE TARGET FROM VARIOUS ARRAY 

CONFIGURATIONS --- 25 

5. FPIB VS. ARRAY APERTURE SIZE - 29 

6. CRITICAL VELOCITY VS. ALTITUDE 32 



7 



LIST OF SYMBOLS 



a energy containing turbulence wavelength, m 

AA incremental area element in the target plane, m 2 

A the area within ’’the Bucket", m 2 

B area of an aperture, m 2 

Cj constant of proportionality, watts m~ 2 

C 2 constant of proportionality, m 0,2 

C L lift coefficient 

d interaperture spacing, m 

D g single turret aperture diameter, m 

D a individual aperture diameter in the array, m 
F array factor 

h altitude, km 

H scale height, km 

I irradiance on the target, watts m ~ 2 

I irradiance on axis on the target for a single aperture, 

watts m -2 

aperture factor 

I 2 radiant emittance of the aperture, watts m -2 

Ij. irradiance at a point in the target plane at grid 

J coordinates i, j, watts m -2 

Jj Bessel function of first order 

k wavenumber, m ~ 1 

L aerodynamic lift force, lbf; boundary layer thickness, m 



S 



l a 

M 

N 

P P 

R 

R 

S 

T 

T o 

V 

X 

a 

Y 
A 2 

0 

K 

X 

A 

v 

y 

p 

Ap 



maximum array dimensions, m 

the number of grid points within "the Bucket"; Mach 
number 

the number of apertures in an array 
power projected, watts 

perpendicular distance between aperture and target 
planes, m 

Reynolds number 

wing surface area, ft 2 

temperature, K 

stagnation temperature, K 

flight velocity, m sec " 1 and knots 

distance from leading edge of flat plate, m 

extinction coefficient for coherent beam in the 
turbulent boundary layer, m ~ 1 

ratio of heat capacities 

geometric path difference between a distance extending 
from a point in the target plane to a reference 
position in the j t * 1 aperture, and the distance R, m 

angular divergence from beam axis 

proportionality constant which relates refractive 
index to air density 

radiation wavelength, m 

turbulence integral scale, m 

kinematic viscosity, m 2 sec 

dynamic viscosity, kg m^sec " 1 

atmospheric density, kg m " 3 

local change in density, kg m ” 3 
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I. INTRODUCTION 



Propagation of laser energy from an aircraft turret 
designed with today's technology to an airborne target is 
unacceptably degraded when the target is in the quadrant 
between 120° - 240° relative bearing due to the turbulent 
boundary layer and separated flow on the aft section of the 
turret. For Surface to Air Missile (SAM) defense, this is 
an especially interesting problem since the most dangerous 
threat to an aircraft is likely to approach from behind. 
Figure la is a diagram of such a turret. The fairing is 
needed to curtail heavy oscillatory aerodynamic loading. 

Were the flow around such a turret laminar, an undis- 
persed beam could be projected into the volume which the 
separated flow presently inhibits. Reduction of turret scale 
could - if sufficient - result in laminar flow. However, at 
aircraft flight velocities and turret scales of practical 
size, a separated flow region will exist in any case 
(Fig. lb). 

A future development may permit undistorted beam propa- 
gation from a window on the aft section of a streamlined 
body. In this event, a similar, though less severe beam 
degradation problem would probably still exist on the large 
scale due to the turbulent boundary layer (Fig. lc) . Note 
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the large scale turret would have a turbulent wake. A laser 
beam propagated through the turbulent wake would be severly 
degraded. 

The small scale streamlined body however, could achieve 
laminar flow (Fig. Id). Such a body with a future technology 
window could project a beam without degradation from separated 
regions or turbulent boundary layers. Beam dispersion in a 
laminar wake is not significant. Reduction of body size 
would however require more apertures to project an amount of 
power equivalent to that from a single larger aperture. An 
array of many such small turrets in laminar flows with total 
aperture area equal to the area of the single aperture may 
be distributed around the aircraft away from aircraft turbulent 
wakes. If properly configured, such an array will achieve 
greater power on a target aft of the aircraft than any 
other of today's technology turret designs, and perhaps 
greater power on a target than the tomorrow's technology 
single- turret design. Henceforth, discussion is limited to 
laser turrets designed using this future "window" technology. 
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small scale large scale I small scale large scale 



TODAY'S TECHNOLOGY (blunt base bodies on fuselage) 





FUTURE TECHNOLOGY (streamlined bodies isolated 

from fuselage) 





Fig. 1: PRESENT AND FUTURE TECHNOLOGY TURRET DESIGNS 
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II. METHOD 



One seeks a quantitative comparison between the per- 
formance of the single larger aperture projecting a beam 
through a turbulent medium and an array of small apertures 
projecting beams through laminar flows which converge on a 
distant airborne target. For identification, the _Single 
Aperture projecting through a Turbulent boundary layer will 
be SAT; the Array of apertures projecting through Laminar 
flows will be AL. 

Isolation of the effects of the turbulent boundary layer 
on the single aperture and of the performance of the array 
in laminar flows is required to determine if greater power 
on the target can be achieved by the array. To do so requires 
the following assumptions: 

1. No power losses in the atmosphere. 

2. No thermal blooming effects. 

3. No jitter. 

4. A perfect pointer/tracker . 

A reference area is defined as that area contained by 
the first zero of irradiance of the Airy disk which would be 
produced on the target in vacuo by the single aperture. 

This reference area is chosen believing that a single aper- 
ture design delivers sufficient destructive energy to a 
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target when projecting through a laminar flow and that energy 
outside of this area will not contribute to target destruc- 
tion or will be insufficient to disable it. This reference 
area will be called "the Bucket". 

Power i_n the Bucket (PIB) for the SAT and the AL will be 
compared. If the PIB for the array is found to be greater 
than that for the SAT, then the array would perhaps be more 
effective against targets approaching an aircraft from behind. 
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III. DIFFRACTION OF LASER BEAMS 



The distribution of monochromatic (laser) irradiance 
on a target is - among other things - the result of diffraction. 

A. SINGLE APERTURE 

For the single circular aperture in the Fraunhofer region 
without any dispersive atmospheric effects the resultant 
distribution of irradiance is the well known Airy disk. 
Calculations show that 83.8% of the incident power arrives 
within the area contained by the first zero of the Airy 
distribution [Ref. 1]. 

B. ARRAY OF APERTURES 

The array possesses maximum dimensions which are larger 
than the single aperture. Thus, the Fraunhofer condition 
may not be satisfied for the array at a target distance 
where it is satisfied for the single large aperture. 

Figure 2 shows an array of apertures where 0^ is the 
source plane, 0 is the aperture plane, and 0£ is the target 
plane . 
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Fig. 2: COORDINATE SYSTEM FOR AN ARRAY DIFRACTION PROBLEM. 

REPRODUCED FROM STONE. [REF. 2] 

For a plane wave incident on the aperture plane the 
Fraunhofer condition is satisfied when the maximum dimension 
of the array creates a phase difference less than 0.1A in 
the target plane. This occurs when either of the following 
conditions are met. 

1. The beam is collimated, and Lf/RA <<1. 

2. The array is focused at the target plane. 

A lens has the effect of bringing the Fraunhofer region 
closer to the aperture plane. It does this by lengthening 
the geometric path length of light passing through the 
central part of the lens such that the Fraunhofer condition 
is satisfied near the focal point of the lens. Phase retard- 
ation across the aperture plane will have the same effect. 
This retardation can hopefully be accomplished by adaptive 
optics so that targets at close distances will be in the 
Fraunhofer region. 
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For simulation on a computer however, it is much easier 

to chose R>>L 2 / X than- to ' insert phase retardation across the 
A 

apertures . 

C. IRRADIANCE ON A TARGET FROM AN ARRAY OF APERTURES 
The irradiance in the target plane is 



I = Cj^F (1) 

where I x is the aperture factor, F is the array factor, and 
Cj is a constant of proportionality. Ij is a function of 
aperture size and shape, distance to and location on the 
target, and the wavelength. For identical circular apertures 



where 

a = (3) 



2 J 2 (ott) 

O 7T 



( 2 ) 



F, the array factor, is the product [Ref. 3] of the 
conjugates of the sum of the complex phases at each point. 
The phases involved are a function of the array geometry, 
distance to and location on the target, and the wavelength. 

N ikA (j) 



F = 




2 



(4) 



J 



where 

^(j) _ -xsin 0 - ysin d> 

2 2 1 2 



(5) 
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Refer to Figure 2 for the geometry and definition of 0 and 

<j) . C is chosen so that the irradiance on axis in the 
2 1 

target plane is N 2 and has the units of watts/m 2 . Substi- 
tuting (2) and (4) into (1) , the irradiance becomes 
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N ikA (j) 


I = c 
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2 J, (air) 

01 T 
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j 



To visualize the distribution of irradiance at the target 
and to calculate PIB, a grid of points (X 2i , Y 2 ^ ) is created 
in the target plane (see Fig. 2). At each point the phases 
and the argument of J i are computed. This yields I at 
that point; I is stored in an array. Those grid points 
within "the Bucket" contribute to the summation of PIB in 
the following manner: 

The PIB for a continuous intensity distribution would be 



PIB 



AL 



= ff I (x,y) dxdy 
bucket 

This is approximated numerically as 

bucket radius 



PIB 



AL 



which becomes 



I.. AA 



1 j 



bucket radius 

PIB = § VV I.. 

AL M / y / V 13 

i j 



(7) 



(3) 



(9) 



20 



if the following further approximation is made: 



AA 



A 



M 



( 10 ) 



This approximation of AA improves as the density of points 
in the grid becomes greater. Calculation of PIB for the 
Airy disk by this method with about 80,000 points in the 
bucket yields 83.5% of the total Power Projected (P^) in 
the bucket. The published value is 83.81. The program is 
in Appendix A. 

The PIB obtained from Equation (9) is compared to the 
total power projected by the array to find the Fractional 
Power in the Bucket (FPIB) . FPIB is the true measure of 
the "goodness" of a system of one or more apertures. 

pib al 

fpib al ■ ^ 



Power projected is found by use of the relation for the on 
axis intensity from a single circular aperture. 

i- -i 



I 



o 



I 

2 



la. 

AR 



( 12 ) 



where I 2 is the intensity of the field in the aperture, B 
is the aperture area, A is the wavelength, and R is the 
distance to the target from the aperture(s). So that for 
an array 



I 



o 



I 

2 




(13) 
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The total power projected by the array is 



P 



P 



X . NB a 



( 14 ) 



Radio astronomers use arrays of antennas; several papers 
[Refs. 4, 5] on radio astronomy were examined for concepts. 
Further, insects have eyes which contain a large number of 
elements in an array. The principals of insect vision [Refs. 

6, 7] were also investigated. 

The irradiance on the target plane was calculated for 
several array configurations; conditions for Fraunhofer 
diffraction existed. Figure 3 shows the FPIB achieved by 
each configuration. The best FPIB achieved for a 10-cm 
aperture was 48% . 

For arrays in the triangular configuration (probable 
candidate for use on aircraft) with array aperture area equal 
to the area of the single aperture, the ratio D./D„ will 
also determine FPIB. Figure 5 shows this trend which is very 
insensitive to interaperture spacing so long as the Fraunhofer 
condition is satisfied; the behavior is true in general for 
all array configurations and is qualitatively apparent from 
Figure 4. 
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Fig. 3: FPIB FOR VARIOUS ARRAY CONFIGURATIONS 




Fig. 3 (continued) 
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Fig. 4a: IRRADIANCE ON THE TARGET FROM A SINGLE 

CIRCULAR APERTURE 



25 





z 



X 




T 




Fig. 4b: IRRADIANCE ON THE TARGET FROM A TRIANGULAR 

3 -ELEMENT ARRAY 
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Fig. 4c: IRRADIANCE ON THE TARGET FROM A CIRCULAR 

9 -ELEMENT ARRAY 
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Fig. 4d : IRRAD IANCE ON THE TARGET FROM A TRIANGULAR 

9 -ELEMENT ARRAY 
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l 7 ig. 5: FPJB VS. ARRAY APERTURE SIZE 



IV. FLUID DYNAMICS RELATED TO TURRETS 



A. SEPARATED FLOW REGIONS 

Behind blunt bodies even at moderate Reynolds numbers 
turbulent separated flow will occur. Thus, the single large 
aperture in a circular turret must project through a region 
of turbulent separated flow such as that shown in Figure la. 

B. WAKES FROM TURBULENT BOUNDARY LAYERS 

Regardless of the body geometry, when the Reynolds number 
is greater than the critical value, a turbulent wake will be 
present behind the body. Propagation through a turbulent 
wake degrades beam quality more than either a separated flow 
or a turbulent boundary layer. 

C. TRANSITION TO TURBULENCE IN THE. BOUNDARY LAYER 

An array of small scale streamlined turrets can maintain 
laminar flow over the many smaller turrets if the geometric 
scale is small enough. An upper bound on the irradiance 
which mirrors can withstand before distortion and failure 
requires an equal area in the array as the single aperture 
to project equal power. The number of array elements then 
depends on the maximum aperture size that still maintains 
laminar flow since it is determined that arrays of large 
apertures achieve greater FPIB than arrays of small apertures. 



30 



The array apertures are located inside of aerodynamically 
streamlined bodies with X A /D A =6 * (see Fig. Id). This design 
is proposed because of transition at a high Reynolds number, 
somewhere above R =2xl0 6 [Ref. 8]. 

C 



V c = 



R v 
c 



R v 
= _£ 



6 D. 



(15) 



Use of equation (15) and typical values of dynamic vis- 
cosity in the atmosphere as a function of altitude [Table 1, 
Ref. 9] yields a critical velocity below which there is 
laminar flow over the turret in the array. 



TABLE 1 



KINEMATIC VISCOSITY VS. 

height (km) 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 



ALTITUDE IN THE ATMOSPHERE 

yxlO 5 (m 2 /sec) 
1.4638 
1.5841 
1.7171 
1.8648 
2.0289 
2.2119 
2.4164 
2.6456 
2.9032 
3.1936 
3.5221 
3.8948 
4.5601 
5.3389 
6.2508 
7.3185 



*For simplicity in these theoretical analyses it is assumed 
that the streamlined turret diameter is equal to the aperture 
(mirror) diameter in Figures lc and Id. 
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ALTITUDE (KM) 



For various aperture geometric scales, Figure 6 shows 



critical velocity at altitudes from sea level to 15 km. 




Fig. 6: CRITICAL VELOCITY VS. ALTITUDE 



Figures 5 and 6 reveal that aperture geometric scale, the 
altitude where laminar flow is maintained, and FPIB for the 
array are all closely linked. It should be noted that an 
aircraft capable of transporting the laser system may not 
be able to fly at the velocity which corresponds to certain 
altitudes and aperture geometric scales. For example, if 
35% FPIB is desired from a triangular array, this requires 
10-cm apertures in the array if D s = 30 cm. However, at an 
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altitude of 5 km transition from laminar to turbulent flow 
occurs at about 145 knots for a 10-cm turret. A large-body 
jet aircraft probably cannot fly at this speed at an altitude 
of 5 km . 

Thus, an estimate of regions of aircraft operability is 
needed; recalling the expression for aerodynamic lift: 

L = %pV 2 SC L (16) 



Then 




2L 



(17) 



is the minimum velocity for an aircraft in a no-stall condi- 
tion. For most large-body jet aircraft L/S is typically 
about 100 lbs per square foot, and a practical maximum for 

C. is assumed to be 1.2. V as a function of altitude 
L m 

been computed and is also shown in Figure 6. 
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V. INFLUENCE OF TURBULENT BOUNDARY LAYER 
ON POWER IN THE BUCKET 



A. SELECTION OF MODEL 

An estimate of PIB for the large single aperture is sought. 
The literature contains theoretical discussion of beam quality 
for propagation through turbulent boundary layers of uniform 
thickness [Ref. 10]. The angular dependence of beam intensity 
for the SAT is estimated by such a model. 

B. ESTIMATION OF PIB FOR SAT 

Sutton has calculated the time averaged angular dependence 
of beam irradiance for several cases of beam diameters (D ) , 
energy containing turbulence wavelength (a) and boundary 
layer thickness (L) for the above model. 

Data points were taken from Figure 2 of Sutton [Ref. 11] 
which plots I/I Q vs. 0D s /X. The trapezoidal rule was used to 
numerically integrate 



PIB, 



SAT 




- ^ 2irr dr 



0 



(18) 



where 



r 




(19) 
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A number representing PIB is obtained for each case of D s /a 
and aL.* The PIB obtained in this manner is converted to 
FPIB by comparison to the result of the above integration 
using I(r)/I Q for the Airy distribution, which is known. 



fpib sat 



PIB 



SAT 



PIB 



AIRY 



8 3.8% 



( 20 ) 



Table 2 shows these results. The program is in Appendix B. 



TABLE 2 

FPIB FOR SELECTED VALUES OF D s /a AND aL 



aL 

1 2 4 10 




*a is nearly equal to A, the turbulence integral scale 
[Ref. 12] 
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To find the region of Table 2 which most closely corres- 
ponds to the case of the single aperture streamlined turret, 
the following equations are used to estimate aL and D s /a. 

The equation for the extinction coefficient in the turbulent 
boundary layer is obtained from Sutton [Ref. 13]. 

a = 2tf<An 2 >A (21) 

where k is the wavenumber, An the change in refractive index, 
and A the turbulence integral scale which Sutton [Ref. 14] 
estimates as 

A = 0.1 L (22) 

L is the path length through the turbulent boundary layer or 
the boundary layer thickness. This length is estimated using 
the equation for the flow of an incompressible fluid in a 
turbulent boundary layer [Ref. 15] beginning at the front 
edge of a flat plate. 

6.2 

»rl (23) 

pV _ 

The variables p and p have the usual meanings, V is the air- 
craft speed, and X is the distance from the leading edge 
(front of the streamlined body in this case) . The beam 
emerges from the rear of the body at a distance perhaps 4/5 
of the body length back from the front, or SD s . A reasonable 
estimate for L is therefore 



L = 0.37 
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L 



0.37 



5 4 y 
P V 



0.2 



0*8 



D 

s 



(24) 



Let C 2 = 0.37 (625u/pV) , which is clearly dependent on 

altitude and velocity. Substitution shows that 

L = C 2 (h , V) D s °’ 8 (25) 



It is known that 



An = K(P-) f 2 

^ oo 



or 



An = K 



SL 



-h/H . 
e £p. 



(26) 



(27) 



where K is the constant in the equation relating refractive 
index to density, i.e., n=l + KQ/P m , H is the scale height, 

Ap is the change in local density, and is the ambient air 
density. From Sutton [Ref. 16] one obtains the relation for the 
last term in Equation (27) . 



Ap 



= 0.1 



1 - 



Pw 



(28) 



where p w is the density of air on the surface of the turret. 
Assuming an ideal gas and constant pressure across the boun- 
dary layer 



Ap 

Poo 




(29) 
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Liepmann and Roshko [Ref. 17] provide an expression for the 
stagnation temperature (T Q = T ) as a function of the ambient 
temperature (Too) , M the Mach number, and y the ratio of 
heat capacities. 




Since stagnation temperature is equated to wall temperature. 
Equation (30) is valid only for Prandtl number of unity. 

From the above equations it can now be shown that 



aL = 2 (4p-) 2 1 0 ' 3 K sl 



-2 h/H 



M 



1 + 



Y-l 






C 2 D 1 / 5 (31) 

2 S 



Typical values for substitution into (31) are H = 7.6 km, 
h = 5 km, D s = 0.3m, y = 1.4, K gL = 5xl0~ 4 , M = 0.3, and 
C^(5 km, 250 knts) = 0.16 m 0 * 2 . These values produce an 
extinction number (aL) of about l.lxlO" 3 . From Figure 2 of 
Sutton [Ref. 18] it is evident that the beam is negligibly 
affected by the turbulent boundary layer and therefore 
very close to the diffraction- limited case regardless of 
the ratio D s /a. 
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VI. CONCLUSIONS 



Laser turrets designed with today's technology may be 
ineffective against an airborne target aft of the aircraft. 

An array of apertures using a future technology window and 
proper phase retardation can achieve reasonable Fractional 
Power in the Bucket (FPIB) on a target approaching an aircraft 
from behind. This FPIB is largely dependent on the aperture 
diameter of the array elements and is relatively insensitive 
to interaperture spacing as long as the Fraunhofer condition 
is satisfied. 

However, possession of the technology for a small scale 
window (see Fig. Id) for use in the array would probably also 
yield the technology for the large window (see Fig. lc) soon 
thereafter. Given the flexibility of both large and small 
window technologies, the single turret shown in Figure lc 
is by far the better. 

Even though the array achieves laminar flow and small 
windows may be easier to develop, such a system incurs the 
price of a collosal fire control problem, phase control 
problem, and the operational restriction of flight velocities 
near the aircraft stall speed. Surprisingly, the single 
turret design is negligibly affected by the turbulent boundary 
layer and therefore nearly diffraction limited upon exit from 
the boundary layer. Though it does create a wake which is 
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turbulent, reasonable advance detection of targets approaching 
from behind should avoid the necessity to try to project 
directly aft through the wake. 
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APPENDIX A (continued) 
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